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Lya versus X-ray heating in the high-z IGM 
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ABSTRACT 

In this paper we examine the effect of X-ray and Lya photons on the intergalac- 
tic medium temperature. We calculate the photon production from a population of 
stars and micro-quasars in a set of cosmological hydrodynamic simulations which 
self-consistently follow the dark matter dynamics, radiative processes as well as star 
formation, black hole growth and associated feedback processes. We find that, (i) IGM 
heating is always dominated by X-rays unless the Lya photon contribution from stars 
in objects with mass M < 10 8 M Q becomes significantly enhanced with respect to the 
X-ray contribution from BHs in the same halo (which we do not directly model), (it) 
Without overproducing the unresolved X-ray background, the gas temperature be- 
comes larger than the CMB temperature, and thus an associated 21 cm signal should 
be expected in emission, at z 11.5. We discuss how in such a scenario the transi- 
tion redshift between a 21 cm signal in absorption and in emission could be used to 
constraint BHs accretion and associated feedback processes. 
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1 INTRODUCTION 

While the investigation of the very high-2 universe, z ~ 
1000, is possible through the Cosmic Microwave Background 
(CMB) radiation and detection of a handful of objects 
with redshift as high as ~ 7 has recently been possible, 
there is a lack of observational data in the redshift inter- 
val 7 ,$ z ,$ 1000. To detect radiation from this interval, 
telescopes with exceptional sensitivity in the IR and radio 
bands are needed. JWSlQ (James Webb Space Telescope), 
for example, with its njy sensitivity in the 1 — 10 ^m in- 
frared regime, is ideally suited for probing optical-UV emis- 
sion from sources at z > 10. Similarly, the planned genera- 
tion of radio telescopes as SK^EI (Square Kilometer Array), 
LOFAI0 (LOw Frequency ARray), 21cmA (21cm Array) and 
MWiQ (Murchison Widefield Array) will open a new observa- 
tional window on the high redshift universe. In particular, 
the detection of the 21 cm line associated with the hyper- 
fine transition of the ground state of neutral hydrogen, holds 
the promise to shed light on the reionization process and its 
sources. 

This observational progress has been accompanied by 
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a flourishing of theoretical activity, aimed at providing pre- 
dictions for the pl anned gene ration of telescopes. It has long 
been known (e.g. iFieldl Il959h that neutral hydrogen in the 
intergalactic medium (IGM) and gravitationally collapsed 
systems may be directly detectable in emission or absorption 
against the CMB at the frequency corresponding to the red- 
shifted HI 21 cm line (associated with the spin- flip transition 
from the triplet to the singlet ground state). iMadau et al.l 
1 19971 ) first showed that 21 cm tomography could provide 
a direct probe of the era of cosmological reionization and 
reheating. In general, 21 cm spectral features will display 
angular structure as well as structure in redshift space due 
to inhomogeneities in the gas density field, hydrogen ion- 
ized fraction, and spin temperature. Several different signa- 
tures have been investigated in the recent literature among 
which fluctuations in the 21 cm line emission induced by 
the "cosmic web" (|Tozzi et alj |2000|). by the neutral hy- 
drogen surviving reionization (e.g. Ciardi fc Madaul 120031 : 



rg 

iFurlanetto eT al. 2004; Mellema et aflfe ooe) and a global fea- 
ture ( "reionization step" ) in the continuum spectrum of the 
radio sky that may mark the abru pt overlapping pha se of 
individual intergalactic HII regions dShaver et al.lll999l ). 

A key feature for the observation of the line in 
emission is that the IGM should be heated above the 
CMB temperature. While pre-heating can happen due 
to, e.g., dark matt e r particles decay o r annihilation (e.g. 
iMapelli et~ai] l200fj ; IValdes et ail 120071 ), the main source 
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of heating are Lya and X-ray photons. Although the im- 
pact of t hese sources has been investigated by several au- 
thors fe g.lMadau et a0l997l:IChen fc Miralda-Escudel | 2004 | : 
Nusseij|2005l; iKuhlen et al.ll2006l; I Ciardi fc Salvaterrall2007l ; 
Pelupessv et al.l 120071 ; iRipamonti et al.1 120081 ). none of the 
above studies has approached the problem considering both 
sources within a self-consistent model. Here, we estimate the 
Lya and X-ray photon production using the simulations by 
I Pelupessv et al.l (|2007T ). which were designed to follow the 
evolution of both quasar and stellar type sources, includ- 
ing their associated feedback effects. The above estimates 
are then used to calculate with a semi-analytic approach 
the evolution of the temperature of an IGM at the mean 
density. 

The paper is structured as follows. In Section [5] we 
describe the main ch aracteristics of the simulations by 
I Pelupessv et al.l (|2007l ). while in Section [3] we give the pre- 
scription used to derive the Lya and X-ray photon count. In 
Section [4] we estimate the IGM heating due to the emitted 
photons, in Section [5] the consequences for the observability 
of the 21 cm line from neutral hydrogen and in Section[S]we 
give our conclusions. 



2 SIMULATIONS OF HALO COLLAPSE AND 
BLACK HOLES GROWTH 

As mentioned in the Introduction, in this paper we mak e 
use of the simulations described in I Pelupessv et all (|2007l ). 
Here we briefly outline the main characteristics and we refer 
the reader to the original paper for details. 

The simulations are designed to investigate the phys- 
ical conditions for the growth of intermediate mass seed 
black holes (BHs; possibly the remnants of a first gen- 
eration of massive stars), using the parallel cosmologi- 
cal TreeP M-Smooth Part icle Hydrodynamics (SPH) code 
GADGET2 (|Springelll2005l ) in the standard ACDM modefl 
The initial conditions for the simulation correspond to 
isolated spherical overdensities ('top-hat') endowed with 
an appropriate Zeldo vich power spectrum (similar to 
IBromm fc Larson! liool ). During collapse of the parent halo 
(all halos considered have T V1I > 10 4 K) the seed holes 
are incorporated through mergers into larger systems and 
accrete mass from the surrounding gas. The interstellar 
medium (ISM), star formation and supernovae feedback as 
well as black hole accretion and associated feedback are 
treated self-consistently by means of sub-resolutions mod- 
els. In particular, the m ultiphase model for star form ing gas 
has been developed bv ISpringel fc Hernauisd (|2003T j . while 
the prescription for accretion and feedback from massive 
bla ck holes is the one de veloped bv lDi Matteo et al.l (120051 ) 
and lSpringel et al. (2005). Technically, black holes are repre- 
sented by collisionless particles that grow in mass by accret- 
ing gas from their environments. The (unresolved) accretion 
onto the black hole is related to the large scale (resolved) 
gas dis tribution using a Bondi-Hoyle-Lyt tleton parameteri- 
zation |Bondilll952i ; lBondi fc Hwlelll944h . and it is limited 
by the Eddington rate. A fraction ei of the radiative en- 
ergy released by the accreted material is assumed to couple 

5 n A = 0.7, n m = 0.3, tt b = 0.04 and H = lOO/i km s" 1 Mpc 
(h = 0.7), where the symbols have the usual meaning. 



thermally to nearby gas. In addition, two black hole par- 
ticles are assumed to merge if they come within the spa- 
tial resolution of the simulation (i.e. within the local SPH 
smoothing length) and their relative speed lies below the 
local sound speed. The seeding procedure constists in se- 
lecting objects with a mass of 10 6 Mq and place a seed BH 
of Msccd = 10 3-4 Mq in them if they do not already contain 
a BH. 

Different collapse scenarios have been investigated by 
choosing the total mass M of the host halos, the redshift 
of collapse z v ir and the spin parameter A. The halos consid- 
ered have a mass of M = 10 8 , 10 9 , 10 10 and 10 11 M Q , with 
a virialization redshift of z v ir = 16, 12, 10 and 7.5, respec- 
tively. A standard spin parameter of A = 0.03 is adopted, 
but simulations with A = or A = .05 have determined n o 
significant difference in the results l|Pelupessv et al.ll2007r i. 
The simulations used here are the highest resolution ones, 
having N = 10 7 particles (in dark matter and gas), with a 
corresponding spatial resolution of a about hundred parsec. 
Two extreme values for the BH feedback have been consid- 
ered: ef = corresponding to the case without any black hole 
feedback (this can be used for comparison) and et = 0.5, the 
value that has been used to reproduce the observed normal- 
izatio n of the local Mbh-Q" re lation in galaxy merger simula- 
tions (|Di Matteo et al1l2005h as well as in full cosmological 
hydrody namical simulations w here the same modeling was 
applied (|Di Matteo et al. 2008). All simulations here are run 
to a final redshift of z ~ 6. 

Despite the analysis of a quite broad parameter space, 
the simulations do not account for possible effects associated 
with gravitational recoil (for a thourough discussion on the 
implications we refer the reader to the original paper). 

The simulations are used to derive the X-ray and Lya 
photon production as described in the following Section. 



3 PHOTON COUNT 

In this Section we discuss our estimate of the X-ray and Lya 
photon production and associated background. Our refer- 
ence run has M scc d = 10 4 Mq, A = 0.03, tt = 0.5. 

The bolometric luminosity associated to a BH accret- 
ing at a rate Mbh{z, Mh) hosted in a halo of mass Mh at 
redshift z is L(z,Mh) = ec 2 Mbh(z, Mh), where e = 0.1 is 
the radiative efficiency for accretion (for a Shakura-Sunyaev, 
and non-spinning black hole) and c is the speed of light. The 
comoving specific emissivity at redshift z is then: 

j BH (u, z) = / r,,',. A z > M h )n PS (z, M h )dM h , (1) 

JM himin \lw{v',z)dv' 

where nps(z, Mh) is the weight for the halo of mass Mh at 
redshift z as computed by the Press-Schechter formalisnjf], 
and Mh,min and Mh,max are the minimum and maximum 
masses for our simulated halo mass, respectively. Mh, m in is 
set by the physics appropriately captured in the simulations, 
which do not follow the chemistry of molecular hydrogen. 
l(v, z) is the ave rage spectrum energy distribution for AGNs 
as computed bv ISazonov et al.l l|2004h . 

6 Before the redshift of virialization of the halo z v i r (Mh), the 
weight has been computed at z„; r . 
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The total cosmic star formation rate density at redshift 



P*(z)= I M±(z,M h )nps(z,M h ), 

>M himin 



(2) 



where M*{z, Mh) is the star formation rate for the halo of 
mass Mh at redshift z. 

For stars, the comoving specific emissivity is computed 

by: 



(3) 



where h(v,t zz i) is the template specific luminosity for a 
stellar population of age t z>z i ( time elapsed betwe e n red - 
shift z' and z) as computed by iBruzual fe Chariot! (|2003l ) 
for Pop II stars wit h Z = 0.2 Z Q and a Salpeter IMF. In 
ICiardi fc Salvaterral (|2007l ) it was shown that, although the 
ionizing photon production is much reduced (by a factor of 
~ 4) compared to metal-free stars with the same IMF, the 
Lya photon emission is similar. For this reason we limit our 
discussion to metal enriched stars. 

The background intensity J(vo, zq) seen at a frequency 
vq by an observer at redshift zq is then given by: 



J{uo, z ) 



(1 + zpf 

4% 



77 

j(u,z)e- T(vo ^^dz, (4) 



where v = vo(l + z)/(l + zo), dl/dz is the proper line ele- 
ment, and j(y, z) — j*(y, z)+jBH(^, z) is the total emissivity. 
t(vo, zo, z) is the optical depth of the medium. 

For Lya phot o ns, r as been computed as in 
ISalvaterra fc Ferraral l|2003l . see Sect. 2.2 for a full descrip- 
tion of the IGM modeling) . The evolution of the Lya back- 
ground is shown in Figure [1] We no te that its intensity 
is mu ch smaller than what found in Cia rdi fc Salvaterral 
l]2007l ). This effect can be attributed to the larger M^^min 
employed here, as set by the smallest simulated halos. 

In the X-rays, 



. . r di 

t{v,Zo,Z) = / —a 

J ZQ 



1 



1+Z0 



ub(z'), 



(5) 



is the photon-ionization cross section per 
a cosmological mixture of H and He 



where u(v) 
baryon of 

l|Zdziarski fc SvenssonllT989l : see also iRipamonti et aL I l2008l 
for a more thorough discussion) and tib(z) = ns(0)(l + z) 3 
is the cosmological baryon number density at redshift z 
(n s (0) ~ 2.5 x 10" 7 cm" 3 ; Spergel et al. 2007). Here, we 
neglect the stellar contribution to the X-ray background 
which is negligible with respect to the AGN one (see also 
lOhlteOOlh . so that j(v,z) = j'bh(^, z). The evolution of the 
X-ray background at v§ — 150 eV (dashed line) and 1 keV 
(dashed-dotted) is shown in Figure [T] Even though a de- 
tailed comparison is beyond the scope of this paper, we note 
that our results here a re overall consistent with those of 
IRipamonti et all (|2008l ) who used a set of semi-analitical 
models to study X-ray heating from early black holes. 

A strong upper limit to the AGN activity at high red- 
shift is set by the u nresolved fraction of t he observed cosmic 
X-ray background jSalvaterra et al]|2007T ). To make sure our 
mini-quasar population does not conflict with this require- 
ment, we compute the contribution of our simulated BHs 
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Figure 1. Evolution of the Lya (dotted line), X-ray at 150 eV 
(dashed) and 1 keV (dashed-dotted) background radiation. The 
solid line is the minimum intensity required for the Lya back- 
ground to be effective in decoupling the spin temperature from 
TcMB- 



to the observed diffuse cosmic X-ray background in the 0.5- 
2 keV and 2-8 keV bands. We find that this contribution is 
~ 0.9 x 10~ 12 and ~ 1.8 x 10~ 12 erg s -1 cm -2 deg -2 , respec- 
tively, that is ~ 4 0% of the most recent estimates of the un- 
resolv ed fraction dHickox fc Markevitcbll2007l : iMoretti et al.l 
2009). Even considering the expected contri bution of faint, 
z < 4 AGNs missed by current deep surveys (|Volonteri et al.l 
2006), our calculation implies values well below the obser- 
vational limits. 



4 IGM HEATING 

Results from previous studies (e.g. Madau et al.ll2004h have 
shown that photons with energies hv > 150 eV typically 
have a mean free path larger than the average separation be- 
tween sources; it is possible therefore to divide the radiation 
field in a fluctuating component with 13.6 eV< hv < 150 eV 
(which creates expanding patchy HII regions), and a nearly 
uniform soft X-ray background component with 150 eV< 
hv < 2 keV (which ionizes the IGM homogeneously). Here 
we concentrate on the effect of the latter component on the 
thermal evolution of the IGM. 

We can calculate the energy rate per H atom (in units 
of [erg s -1 ]) at redshift z as: 



<.2kcV 

E x {z) =4vr / J{E, 

Jl50cV 



)a(E)dE, 



(6) 



where J(E, z), with E — hv, is computed by eq. [4] 

Of the above energy, a fraction fheat goes into heat- 
ing of the IGM and a fraction fi on goes into H ionization 
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(|Shull fc van Steenberdfl985l ; IValdes fc Ferrarall200Sl ) 0. We 
have verified that by using the two parameterizations of 
fheat the overall discrepancy in our results (in terms of IGM 
temperature) is at most a few percent . In th e following, we 
adopt the results bv lValdes fc Ferraral I^OOct ) . 

The evolutio n of the IGM tem p eratu re, Tigm, is 
followed as in ICiardi fc Salvaterral (120071 ) (based on 
IChuzhov fc Shapirol l2007h. i.e. Tigm is calculated with 
RECFAST ( Seager et al.lll999l ) until the first sources of radi- 
ation turn on. Subsequentely, taking into account the pres- 
ence of the heating from X-ray photons and neglecting the 
contribution from deuterium, which turned out to be negli- 
gible, the evolution of Tigm is regulated by: 



dTiGA 



3k 



(H a + H x ) 



4Tigm 

3t ' 



(7) 



Here H a and H x — fheat{z, xun)E x (z) are the heating rates 
per H atom due to Lya and X-ray phot ons, respectively. We 
refer the reader to the original paper (ICiardi fc Salvaterral 
120071 ) for the derivation of H a . N a (the number of photons 
that pass through the Lya resonance per H atom per unit 
time) necessary to calculate H a has been derived from the 
model described in the previous Section. 

It should be noted that the impact of Lya pho- 
ton scattering on the evolution o f the IGM temperature 
has b ee n recen t ly rev i sed by e.g. IChen fc Miralda-Escudel 
ll2004li. Iffiratal J200rj|. IPritchard fc Furlanettd (120061) and 
IChuzhov fc Shapirol l|2007h . lChen fc Miralda-Escudel (|2004f ) 
included atomic thermal motion in their calculations, find- 
ing a heating rate se veral orders of magni tude lower than the 
previous estimate bv lMadau et~ai1 (|l997f ) . In addition, while 
"continuum" photons (with frequency between the Lya and 
Ly/3) heat the gas, "injected" photons (which cascade into 
the Lya from higher atomic resonances) cool the gas, re- 
sulting in an effective cooling from Lya photons at tem- 
peratures above 10 K. IChuzhov fc Shapiro! l|2007h though, 
showed that the cascade which follows absorption of pho- 
tons in resonances higher than the Lya happens via the 2s 
level rather than the 2p level. Thus, the number of "injected" 
photons and their cooling efficiency are re d uced compared 
to the estimate of IChen fc Miralda-Escudel |2004) and Lya 
photon scattering can be an efficient heating source also at 
tempe ratures higher than 1 K. H ere we follow the calcula- 
tion of IChuzhov fc Shapirol i|2007h . 

As H x is also a function of the ionization fraction, eq.[7] 
needs to be solved together with the following equation: 

— = fion(t,xmi)Nx—^- + k co un e xm — k rec n e xmi, (8) 
at riHi 

where xmi and xm are the HII and HI fraction respectively, 
riH and nm are the number density of H and HI respectively, 
kcoii and fc rec are the collisional and recombination rates, N x 
is the rate of X-ray photons per H atom. 

In Figure [2] the evolution of the IGM temperature as 
given by eqs.[7]and[8]is shown, together with the CMB tem- 
perature, Tcmb, and the IGM temperature, Tk, in the ab- 
sence of any heating mechanisms. The solid curve includes 
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Figure 2. Temperature evolution as a function of redshift: total 
(solid black line), total in the absence of black holes feedback 
(long-dashed green line), total with an extension to lower masses 
of the contribution of stellar type sources (long-dashed-dotted 
cyan line) as determined by X-ray heating from black holes only 
(dotted red line) and by Lya photons from stars only (short- 
dashed-dotted blue line). The upper (lower) dashed curve is the 
CMB temperature (the IGM temperature in the absence of any 
heating mechanisms). 



the contribution to the heating of X-rays from BHs, Lya 
photons from both BHs and stars. It is clear from the Fig- 
ure that the contribution from the Lya photons is irrelevant 
for the evolution of the IGM temperature, which is instead 
dominated by the X-ray heating (see Sec[6]for a more exten- 
sive discussion). We note that the contribution to the Lya 
photons from the BHs is negligible compared to that from 
the stars. 



5 SPIN AND DIFFERENTIAL BRIGHTNESS 
TEMPERATURE 

As the physics behind the emission/absorption of the 21 cm 
line has been di scussed extensively by several authors (for 
a review see e.g. iFurlanetto et al.ll200rj). here we just writ e 
the relevant equations. following lCiardi fc Salvaterral l|2007h . 
T he evolution of the spin te mperature, T s , can be written 
as IChuzhov fc Shapirol (120061 ): 



Tcmb + (ya, e ff + J/ c )7igm 



1 + Va.eff + 2/c 

where Tcmb is the CMB temperature; 



Vc = 



-(C n + C e + C p ) 



0) 



(10) 



2 ) and fi„n(z,mn) =0.3846(1 



„0.5420U.1952 



^4loTlGM 

is the coupling efficiency due to collisions with H atoms, 
electrons and protons, T* = 0.068 K is the temperature 
corresponding to the energy difference between the singlet 
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and triplet hyperfine levels of the ground state of neutral 
hydrogen and Aio = 2.85 x 10~ 15 s _1 is the spontaneous 



emission rate. For the de-excitation rates Ch, C, 
we have adopted the fits used by Kuhlen, Madau 



and C p 
: Mont- 
gomery (2006; see a lso th e original papers bvlSmithl|l966l; 



lAllison fc Dalgarnol 1 19691 ; lLisztll200l] ; IZygelmanl 120051 ). i.e 



Ch = iihk, C e = n e 7 e and C p = 3.2n p K with riH, n e , n v 
hydrogen, electron and proton number density, « = 3.1 X 
lO- 11 ^ ^ 7 exp(-32/TiGM) cm 3 s" 1 and log^e/lcmV 1 ) = 
-9.607 + 0.51og(ri G M)exp[-(logTi G M) 4 - 5 /1800]. Finally, 



-0.3(1 + 2) 



1/2—-2/3 



1 



0.4 



V 



(11) 



Tigm / 

is the effective coupling efficiency due to Lya scattering 
which takes into account the back-reaction of the reso- 
nance on the Lya spectrum and the effect of resonant pho- 
tons other than Lya, and y a = PioT*/(AioTigm)> with 
Pro ~ 10 9 J a s" 1 radiative de-excitation rate due to Lya 
photons. 

In the upper panel of Figure [3] the spin temperature is 
shown as a function of redshift together with Tomb and the 
IGM temperature in the absence of heating mechanisms Tk 
(upper and lower dashed lines, respectively). As the inten- 
sity required for a Lya background to be effective in decou- 
pling the spin temperature fro m Tqmb is J a ^ 10~ 22 ( 1 + 
z) ergs cm" 2 s" 1 Hz" 1 sr" 1 (e.g. ICiardi fc Madau|[2003l , see 
Fig-HJ, T s gets coupled to Tigm as early as z ~ 19, but only 
at z ~ 11 does it become larger than Tomb due to X-ray 
heating. 

The differential brightness temperature, 8Tb, between 
the CMB and a patch of neutral hydrogen with optical depth 
r and spin temperature T s at redshift z, can be written as: 



ST h = (1 



- T \ Ts — Tea 



' 1 + z 

For a gas at the mean IGM density, the optical depth is 
3Ai A 3 T* mnO) 



(12) 



r(z) 



(13) 



32tt T s H(z) ' 

where A is the wavelength of the transition, nm is the num- 
ber density of neutral hydrogen and H(z) is the Hubble 
constant. 

In the lower panel of Figure [3] 5Tb is plotted for the spin 
temperature calculated above, reflecting its behavior, with 
a signal in emission (absorption) expected for z ^ ( ^ )11. 



6 DISCUSSION AND SUMMARY 

In this paper we have estimated the effect of X-ray and Lya 
photons on the IGM temperature. This was done calculat- 
ing the photon production from a population of both stars 
and quasars based on the self-consistent modeling of these 
so urces and associa t ed fee dback processes in the simulations 
bv lPelupessv et ail (|2007f ) . Such photon production has then 
been used to determine the evolution of the IGM tempera- 
ture with a semi-analytic approach. The calculations assume 
a gas at the mean density whose ionization history is reg- 
ulated by the effect of the X-ray photons. Thus, they are 
strictly valid as long as the filling factor of gas ionized by 
UV photons is small or in regions of the IGM that have not 
been reached by UV-ionizing radiation. For reference, the 
filling factor of ionized regions, under the assumption that 
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Figure 3. Upper panel: Evolution of the spin temperature (solid 
line). Tcmb and the IGM temperature in the absence of heating 
mechanisms, T^, are plotted as upper and lower dashed lines, 
respectively. Lower panel: Evolution of the differential brightness 
temperature. Lines are the same as in the upper panel. 



it is independent from the effect of the X-rays, assuming a 
gas clumping factor of 10 and an escape fraction of ionizing 
photons from stellar type sources of 10% is ~ 0.004, 0.07, 1 
at z = 15, 10, 7.1. It should be noted here that the bulk of the 
UV photons is produced by stars rather than BHs, while the 
contribution to the ionization fraction from X-rays is only a 
few percents. 

The results presented of course depend quite strongly on 
the parameters and assumptions adopted in the simulations 
described in Section [2] As already mentioned, while e.g. the 
value of the spin parameter has a negligible effect on the 
final outcome, the presence of AGN feedback (but not the 
specific value of e/, as long as 7^ 0) is crucial. In our reference 
run we have used a vale ef = 0.5, which corresponds to 
coupli ng only 5% of the radiated energy. In lPelupessv et all 
(2007) we have shown that no significant difference in the 
BH accretion rates are seen for £/ = 0.5, but that of course 
cases with e/ = can imply much larger accretion rates (and 
associated X-ray photons). For the sake of illustration, we 
have performed the same calculations also for the extreme 
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case without feedback, i.e. for e/ = 0. Note that for our 
calculation this is a completely unrealistic model, which we 
consider only as a parametric study. The result is shown 
in Figure [2] as long-dashed line. As expected, the value of 
the IGM temperature raises much more quickly than in the 
reference case and becomes larger than Tcmb already at 
z ~ 16. This would result also in an earlier 21 cm signal in 
emission. This case is both unrealistic and also results in a 
contribution to the unresolved X-ray background which is 
orders of magnitudes larger than the available observational 
limits and thus can be discarded. 

So it remains that our results depend most strongly on 
the choice of Mh,min and Mh, ma x in eqs. Q] and [2] which in 
our case we have taken from the masses of the simulated 
halos. While larger values of Mh, ma x are expected to have 
a negligible impact because of the paucity of such objects 
at the redshifts of interest, a smaller Mh,min (and its asso- 
ciated high abundance) would result in higher photon pro- 
duction. While extending eq.[2]to smaller masses is straight- 
forward, it may not be appropriate as the contribution from 
these low mass halos depends on more complicated forms 
and strength of the feedback active in t he relevant range 
of re dshift (for a review on feedback see ICiardi fc Ferraral 
2005). Although a consensus on the detailed effects of feed- 
back (in particular radiative feedback) on small structure 
formation has not been reached yet, most authors agree 
that the presence of feeback delayes or suppresses the for- 
mation of small mass objects, and thus associated star 
formation, depending on a variety of physical quantities 
such as reds hift, halo mass and molecular hydrogen con- 
tent (see e.g. Machacek et alj 2001 L Susa fc Umemurall2006l ; 



lO'Shea fc Normanl 120071 : IWhalen et all l200Sh . A proper 
treatment of such feedback though is beyond the scope 
of the present paper. Similarly the emission properties of 
mini-QSOs hosted in such small halos is largely unknown 
and may not be described by the average AGN spectrum 
here adopted. Thus, we have neglected the contribution of 
these objects to the total X-ray photon production (see 



Madauetal ] 120041 : iDiikstra et al ] 120041 : ISalvaterra et al.l 
20051 for a discussion about the effect of mini-QSOs in the 



early universe). 

As a reference though, an upper limit for the contribu- 
tion of Lya photon to the IGM temperature is plotted in 
Figure [2] Here (long-dashed-dotted cyan line) we show the 
effects accounting for the contribution to the star formation 
rate of objects with mass as small as 10 7 Mq. As expected, 
in this case, the IGM temperature starts to raise at an ear- 
lier time due to Lya heating and becomes larger than Tcmb 
at z ~ 12, while X-ray photons dominate heating at z 10. 
Note however, this considers that these halos will not pro- 
duce any X-ray photons, i.e. the accretion onto their black 
holes must be completely quenched. This is likely a con- 
flicting requirement, as both star formation and black hole 
accretion will somewhat depend on the same gas supply. To 
summarize this point, in order to have 21 cm line in emis- 
sion, X-ray heating from BHs and/or extremely efficient star 
formation (albeit, as discussed above, somewhat unfeasible) 
in halos with Mu < 10 s Mq is required. The heating due to 
stars only in larger halos as shown in Figure [5] is however 
unable to bring the gas temperature above Tcmb before the 
reionization is completed. 



To summarize the main results of this work, we find 
that: 

• unless the Lya photon contribution from stars in small 
mass objects is somehow strongly enhanced in comparison 
to the X-ray contribution from BHs in the same objects, the 
IGM heating is always dominated by X-rays; 

• in this case, the transition redshift between a 21 cm sig- 
nal in absorption and in emission could be used to constraint 
BHs accretion properties; 

• for a case consistent with the available limits on the 
unresolved fraction of the X-ray background in the 0.5-2 keV 
and 2-8 keV bands, the IGM temperature becomes larger 
than the CMB temperature (and thus an associated 21 cm 
signal should be expected in emission) at z 11.5. 
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